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Since any fermionic operator yj can be written as yb = q+ip, where q and p are hermitian operators, 
we use the eigenvalues of q and p to construct a functional formalism for calculating matrix elements 
that involve fermionic fields. The formalism is similar to that for bosonic fields and does not involve 
Grassmann numbers. This makes possible to perform numerical fermionic lattice computations that 
are much faster than not only other algorithms for fermions, but also algorithms for bosons. 
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Since yb 2 = for a fermionic field operator yb, this 
operator does not possess nonzero c-number eigenvalues. 
Therefore, the functional formalism for fermions involves 
anticommuting Grassmann numbers JIJ . This makes the 
notion of a classical fermion field conceptually difficult. 
This is also problematic from a practical point of view, 
because the functional integrals over Grassmann num- 
bers cannot be performed numerically on a computer, 
which makes lattice calculations with dynamical fermions 
difficult. For typical Hamiltonians, the integration over 
fermionic fields can be reduced to the computation of 
the determinants of huge matrices, which costs a lot of 
computer time. Therefore, most of the existing lattice 
calculations are based on a quenched approximation, in 
which fermions are not treated dynamically (see, e.g., j^] 
and references therein). Various methods to improve and 
speed up the treatment of fermions on a lattice have been 
proposed, such as those in Refs. ||, ||, ||, ||, 0]. (Ref. 
contains a longer list of references for various methods.) 
However, the efficiency of all these existing methods can- 
not be compared with the efficiency of methods for func- 
tional integration over bosonic fields ||. 

In this Letter, we construct an entirely new method 
of functional calculus for fermionic fields. It is based on 
the fact that any linear operator yb can be written in a 
unique way as a sum of a hermitian and an antihermitian 
operator. In other words, it can be written as yb — q + ip, 
where q and p are hermitian operators. Since hermi- 
tian operators can always be diagonalized, we can use 
the eigenvalues and eigenstates of q and p to construct 
a functional formalism similar to that for bosons. We 
shall see that q and p play a role similar to that of the 
canonical coordinate and the canonical momentum for 
bosons. This also suggests that q and p can be naturally 
interpreted as observables for fermionic fields, which may 
be relevant for a deeper understanding of the physical 
meaning of quantum and classical fermionic fields. How- 
ever, we emphasize the practical utility of our method 
because the functional integration over fermionic degrees 
of freedom is not only very similar to that for bosons, 
but the numerical computation is even much faster than 
that for bosons. This is because, for each bosonic degree 
of freedom, the eigenvalues q and p can take any real 



value, so lattice calculations require Monte-Carlo simu- 
lations with a large number of different real values. On 
the other hand, it appears that, for each fermionic degree 
of freedom, q and p can take only two different values. 

First, consider one fermionic degree of freedom de- 
scribed by the operators yb and ft that satisfy the anti- 
commutation relations 



{ftft = {ft,ft} = 0, {ftft} = l. 



(1) 



This algebra can be represented on a 2-dimensional uni- 
tary space spanned by the orthonormal basis vectors |0) 
and |1) that satisfy yb\0) = 0, ft\l) = 0, -0+|O> = |1), and 
■011) = |0). The operators vb and ft can be written as 



yb 



ip, 



ft 



q-ip, 



where 



(ft + ft/2, p = i(ft-ft/2. 



(2) 



(3) 



are hermitian operators. They satisfy the anticommuta- 
tion relations 



{q,q} = {p,p} = 1/2, {q,p} = 0. 



(4) 



In particular, we see that q 



,;'2 



V 



1/4, which implies 



that the only eigenvalues of q and p are q± = ±1/2 and 
p± = ±1/2, respectively. The corresponding eigenvectors 
are 

\q ± ) = (|0) ± |1))/V2, \p ± ) = (|0) ± *|1»/V2, (5) 

with the properties q\q) = q\q), p\p) = p\p), for q (p) 
equal to q + or q_ (p + or p_ ) . These eigenvectors satisfy 
the orthonormality relations 



(q\q') = Sqq>, (pb') = <W> 

and the completeness relations 

5>><?i = i, 5»<p| = i. 



They also satisfy 



-nrpq 



(0) 



(7) 



(8) 
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which is similar to the analogous equation for bosons. 
(In order to obtain even closer similarity to bosons, one 
can remove the number ir in the last expression in (||) by 
an appropriate rescaling of the operators q, p, and their 
eigenvalues.) 

Let us now apply the above equations to a calcula- 
tion of the partition function Z = Tr exp(—/3H), where 
H = ujip^ip is the Hamiltonian and T = 1//3 is the 
temperature. Since ijv^n) — n\n) for n equal to 
or 1, the simplest way to calculate Z is as follows [§: 
z = £„H exp(-/9#)|n> = exp(-/3u;0) + exp(-/3uT) = 
1 + exp(— (3lo). However, this simple method does not 
work for a more general case with many degrees of free- 



dom and complicated interactions. Therefore, we present 
a different method, completely analogous to the func- 
tional method for bosons Q, that can be generalized to 
more general cases. We first write the partition function 

as 

z =5>ie-^> =5>i (rK A ^) i«>. w 

i i y =1 / 

where At = 0/N. The next step is to insert the factor 1 
in the forms of (^) 27V times and to put N to infinity, so 
as to obtain 



N 

7 = J^E I IIEE ] (QlP(N)MN)\e 

1 \j = l q(j) P{j)j 
-At Hi 



At H 



q(N))(q(N)\p(N-l)) 



x<K^V-l)|e- A ^k(7V-l))---<g(2)b(l))(p(l)|e- ATif k(l))< g (l)|g). 
Since the Hamiltonian is H = uiip^ip — w(l/2 — 2ipq) = H(p, q) and since At — ► 0, we can write 

(p(j)\e- ATHm kU)) = (P(m - ArH(p,q)\q(j)) = (p(j)\q(j)) e" Ar H W>'««> . 



(10) 



(11) 



Therefore, using (|1 
as 



and 



e f n 5 ee) W 



101) can be written 



(12) 



where 



iv 



-™ PU) + H(p(j),q(j)) 



i=i 

(13) 

and q(N + I) = q. The summation over q is easily per- 
formed, so the final expression is 



lim 



N- 



n^EE 

7 = 1 P(j), 



(14) 



where the periodic boundary condition q(N + 1) = q(X) 
is understood. Similarity to the bosonic case || is ob- 
vious. As we have already discussed, the factor 7r in 
( |l3| ) can be eliminated by an appropriate rescaling. The 
only important difference is the fact that for bosons 
one has to make the replacement (1/2) J2 q (j) ^2p(j) ~* 
(l/27r) / dq(j) J dp(j). Since for fermions each q(j) and 
p(j) can take only two different values (±1/2), we see 
that the numerical computation (with a large but finite 
N) for fermions is much faster than that for bosons. Note 



also that one can formally write (1 
notation as Z = J J T>q(T)'Dp{r) e~ 





with the continuum 
, where r = (3j/N, 



Vq{T)Vp{T) = f[\Y,^ 



t=0 9(t) P (t) 



dr 



-ITT— p(r) 

dr 



H(p(r),q(T)) 



(15) 



(16) 



but here the derivative dq(r)/dT = (q(j + 1) — q(j))/Ar 
is not well defined because it is always equal to or ±oo 
as At — > 0. However, this is not a problem when S is 
calculated as in (|l3|). 

Let us now generalize the analysis to an arbitrary num- 
ber of fermionic degrees of freedom. We label all degrees 
of freedom by a = 1, . . . , L, where a comprises various 
types of degrees of freedom, such as different space points, 
spinor indices, and color indices. The anticommutation 
relations arc 



{Vv <M = $} = o, $} = s a 



(17) 



We can represent this algebra in terms of direct products 
of representations for each pair ip a , The orthonormal 
basis for the 2 L -dimensional representation space can be 
chosen to be the set of 2 L states of the form 



\n a ) = Tii) 



(18) 
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where each n a is equal to or 1. We take the action of 
the operators on these states to be defined by 



A l ---A L \n 1 ) 



A 1 \n l ) 



■>A L \n L ), (19) 



where A a = A a (if; a , ijj a ) with A a being an arbitrary func- 
tion. To determine the action of a differently ordered 
product of operators, one has to use the anticommuta- 
tion relations ( fl7|) to write the action in terms of opera- 
tors ordered as in (|l9|). For example, i/^i \ni) <8> \n 2 ) — 

—4>\ip2 \ n i) ® Vi>i) = ~V'il n i) ® ^2\ n 2)- The hermitian 
conjugation of (|19l ) determines the action of the operators 
on the left: 

(ni|®- (m| A\---A\ = (n L \$ L ®---®{ni\A\. (20) 

To construct the functional formalism, we write ip a — 
Qa + ipa, so that q a and p a are hermitian operators that 
satisfy 

{Qa,qb} = {Pa,Pb} = (l/2)5 ab , {q a ,p b } = 0. (21) 

We introduce the states \q) = a |g a ) and \p) — a \p a ), 
which have the properties q a \q) = q a \q) and p a \p) — 
p a \p)- These states satisfy the orthonormality and com- 
pleteness relations analogous to (^) and (|7|), respectively, 
while the generalization of (§) is 



(P\q) = 



-nrpq 



(22) 



where pq = J2 a PaQa- 

We want to calculate the partition function for a 
Hamiltonian of the form 



H = ^H {1) 4> + {i>H {2) ^ + h.c), 



(23) 



where $ H^ip = b ipb, and similarly for the 

second term . We can write 



H a b ~ A a b + S a b, S a b — S a b + S a b 



(24) 



where A a b — (H„ b — H^)/2 = —Ab a is the antisymmet- 
ric part, S ab = {H^ + Hj ) 1 a > )/2 = S ba is the symmetric 
part, S a b = S a S a b, and S' aa = 0. Using the anticommu- 
tation relations (pl|), we can write 

J2^ H ab^ = (1/2) X)S."2i 

a,b a a 

+ ^q a A ab qb + ^p a A a bPb- (25) 

a,b a,b 

We want the q a operators to be ordered as in (|l^) , so we 
write 



^q a A ab qb = qaAgbqb + ^ qaA a 



bQb 



^ q a A a bqb + ^ QbAbaqa 
a<b a<b 

2 ^ QaA ab qb 



a<b 



a,b 



where 



Aab, a <b, 
■ h II. a ■ b, 



A; 



0, a < b, 
A ab , a> b. 



(26) 



(27) 



Similarly, we want the p a operators to be ordered as in 
(pp|), so we write 



^PaAabPb = 2^2p a A> b p b . 



(28) 



a , b 



a . b 



We also apply a similar procedure for the second term 
in ((H). Without losing on generality, the matrix is 
taken to be antisymmetric. Defining the matrices A^ 
and A(~^ as A a ^ = ±H^* , we can write the Hamil- 
tonian ( p3| ) in the final form as 

H(p,q) = ^TrS- 2ipSq + 2qA<q + 2pA>p 

+2qA<--^<q - 2pA { ->p + 2ipA^q. (29) 

This form of the Hamiltonian has the property 

(p|£f(p,ff)| g ) = (p|g)fr(p, ? ). (30) 

Therefore, we can repeat the steps that led to ( |l~4] ) to 
obtain 



N 



lim 

N— >oo 



n^EE 

qU) p(J)j 



(31) 



a . b 



a<b 



a>b 



where S has the form (Uj) with the Hamiltonian H(p,q) 
given by (§). In &7E qU) = U a E qaU) , E pU) = 
Ila ^2p a (j) ■ 

Needless to say, the formalism developed above can 
also be applied to the calculation of the partition func- 
tion at finite chemical potential. Moreover, essentially 
the same formalism can be applied to the calculation of 
correlation functions among fields at different times, in a 
way similar to that for bosons. Using the fact that H is 
the generator of infinitesimal time translations, one finds 
similar expressions with At — ► iAt. 

Concerning the application of our formalism to con- 
crete physical models, two additional remarks may be in 
order. First, one should be careful when the Hamilto- 
nian is written in terms of Dirac fields ip , ip ^ . For 
example, if there are two degrees of freedom, one corre- 
sponding to a particle and the other to an antiparticlc, 
then the "Dirac field" is ^jp = + $J)/\/2. Second, 
although the operator N = ^ a i ijr a ip a has the property 



4 



N\n) = (X^a n o)l n )' the operator of the physical num- 
ber of particles, in general, has a more complicated form 
^Vphys = Y, a ,b^l N ab^b- For the explicit form of N ab in 
the continuum case, see JToj ] . 

In this Letter, we have presented a general formalism 
rather than apply it to concrete physical problems. How- 
ever, it is straightforward to apply the developed for- 
malism to any physical model, such as lattice QCD. We 
strongly believe that this formalism will be widely ap- 
plied to lattice QCD and lead to a number of important 
phenomenological results that previously could not be 
obtained because they required too much computer time 
with previous algorithms for fcrmions. 

Finally, let us shortly discuss the representation of 
quantum fcrmionic states in terms of wave functions 
*&{q) written in the ^-representation. Note that it is al- 
ready known how to represent fermionic quantum states 
in terms of wave functions ^(ip), where -0 is a Grass- 
mann variable However, ^(ip) is a Grassmann val- 
ued quantity, so ty*(tp)^(ijj) cannot be interpreted as a 
probability density. This problem does not exist for ^(g). 
We first introduce the base functions 

*+(?) = (q\q+) = *-(?) = (q\Q-) = ©(-<?), 

(32) 

where &(q) is the step function equal to 1 (0) for q = 
1/2 (q = -1/2). From (|) and fll we see that p\q±) = 



anticommutation relations (|j) and ([!]), having in mind 
that the operators act on functions f(q) defined for q G 
{—1/2, 1/2}. For time dependent wave functions, i.e., in 
the Schrodinger picture, the Schrodinger equation is 



H(-iqe q ,q)V(q,t) =id t V(q,t). 



(35) 



The expected value of an arbitrary operator O = 0(q,p) 
is given by 

(6(t)) = ]T **(<?, t)0(q, -iqe q )V(q, t). (36) 



The generalization of Eqs. fl32j)-(|36j) to the case with 
many fermionic degrees of freedom is obvious. 
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